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1. Introduction 

The Chermnykh's problem is new kind of restricted three body problem which was first time 
studied by (Chermnykh 1987). (Papadakis & Kanavos 2007) given numerical exploration 
of Chermnykh's problem, in which the equilibrium points and zero velocity curves studied 
numerically also the non-linear stability for the triangular Lagrangian points are computed 
numerically for the Earth-Moon and Sun-Jupiter mass distribution when the angular velocity 
varies. The mass reduction factor gi = 1 — ^ is expressed in terms of the particle radius a, 

density p radiation pressure efficiency factor x (in C.G.S. system) :gi = 1 — ^"^^ — X- Where 
Fp is solar radiation pressure force which is exactly apposite to the gravitational attraction 
force Fg and change with the distance by the same law it is possible to consider that the 
result of action of this force will lead to reducing the effective mass of the Sun or particle. 

(Ishwar & Kushvah 2006) examined the linear stability of triangular equilibrium points 
in the generalized photogravitational restricted three body problem with Poynting- Robertson 
drag, L4 and L5 points became unstable due to P-R drag which is very remarkable and 
important, where as they are linearly stable in classical problem when < < jJiRouth — 
0.0385201. Further the normalizations of Hamiltonian and nonlinear stability of ^4(5) in the 
present of P-R drag has been studied by (Kushvah et al. 2007a, Kushvah et al. 20076, Kushvah 
et al. 2007c) 

Normal forms are a standard tool in Hamiltonian mechanics to study the dynamics in 
a neighbourhood of invariant objects. Usually, these normal forms are obtained as divergent 
series, but their asymptotic character is what makes them useful. From theoretical point 
of view, they provide nonlinear approximations to the dynamics in a neighbourhood of the 
invariant object, that allows to obtain information about the real solutions of the system by 
taking the normal form upto a suitable finite order. In this case, it is well known that under 
certain(generic) non-resonance conditions, the remainder of this finite normal form turns 
out to be exponentially small with respect to some parameters. Those series are usually 
divergent on open sets, it is still possible in some cases to prove convergence on certain sets 
with empty interior(Cantor-like sets) by replacing the standard linear normal form scheme 
by a quadratic one. 

From a more practical point of view, normal forms can be used as a computational 
method to obtain very accurate approximations to the dynamics in a neighbourhood of the 
selected invariant object, by neglecting the remainder. They have been applied, for example, 
to compute invariant manifolds or invariant tori. To do that, it is necessary to compute the 
explicit expression of the normal form and of the (canonical) transformation that put the 
Hamiltonian into this reduced form. A context where this computational formulation has 
special interest is in some celestial mechanics models, that can be used to approximate the 
dynamics of some real world problems. 

The linearization of Hamiltonian, have been studied numerically and analytically. We 



2 



have examined the effect of gravitational potential from the belt, oblateness effect and 
radiation effect. 



2. Equations of Motion and Position of Equilibrium Points 

Let us consider the model proposed by (Miyamoto & Nagai 1975), according to this model 
the potential of belt is given by: 

b^J4|a.- + (a + 3A')|(a + .Vf 

[r^ + {si + N)Y 

where M5 is the total mass of the belt and = + y'^, a, b arc parameters which determine 
the density profile of the belt, if a = b = then the potential equals to the one by a point 
mass. The parameter a "flatness parameter" and b "core parameter" , where = V^;^ + b^, 
r = a + b, z = 0. Then we obtained 

and = ^;i^|ij372, Vy = ^2!^2f/2 ■ As in (Kushvah 2008a, Kushvah 2008&), we consider 
the barycentric rotating co-ordinate system Oxyz relative to inertial system with angular 
velocity u and common 2;-axis. We have taken line joining the primaries as x-axis. Let 
mi, 7712 be the masses of bigger primary (Sun) and smaller primary (Earth) respectively. Let 
Ox, Oy in the equatorial plane of smaller primary and Oz coinciding with the polar axis of 
m2- Let Ve, Vp be the equatorial and polar radii of m2 respectively, r be the distance between 
primaries. Let infinitesimal mass m be placed at the point P{x,y,0). We take units such 
that sum of the masses and distance between primaries as unity, the unit of time taken such 
that the Gaussian constant of gravitational = 1. Then perturbed mean motion n of the 
primaries is given by = 1 + ^ + tt^^W^' where = (1 - A*)gi^^ + ^2 = 
is oblateness coefficient of mo- where a = — '-^ — is mass parameter, 1 — u = — '-^ — with 

^ " m\+m2 ' ^ mi+m2 

mi > 7712. Then coordinates of mi and m2 are (a;i,0) = (— //,0) and (0:2, 0) = (1 — /u, 0) 
respectively. The mass of Sun mi ^ 1.989 x 10^"^% ^ 332, 946m2(The mass of Earth), hence 
mass parameter for this system is /jl — 3.00348 x 10~^. In the above mentioned reference 
system, we determine the equations of motion of the infinitesimal mass particle in xy-plane. 
Now using (Miyamoto & Nagai 1975) profile and (Kushvah 2008a, Kushvah 20086), the 
equations of motion are given by: 

X - 2ny = C4 - T4 = Q^, (3) 
y + 2nx^Uy-Vy^ Qy (4) 



where 



„ _ 2 {1 - fJ,)qi{x + fj,) _ ijl{x + n - 1) _ 3 ijlA2{x + n - 1) _ M^x 
\ix — n X -3 -3 - -g 



1 



2 r| (^2 + 2^2)3/2' 
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[I - ii)qiy ny 3 nA2y 



Mby 



rl 2 r2 



r2 + T2) 



3/2 ' 



2 ri r2 2r| (^2 _j_ 2^2^1/2 

The energy integral of the problem is given by C = 2Q — x'^—y'^, where the quantity C is the 
Jacobi's constant. The zero velocity curves C = 2Q{x,y) are presented in figure ([I]) for the 
entire range of parameters742 , M;, and qi = 1,0. We have seen that there are closed curves 
around the L4(5) so they are stable but the stability range reduced( [2])] due radiation effect. 
The closed curves around L^^^ disappeared when gi = 0. 



(5) 



The position of equilibrium points are given by putting = fly = i.e., 



ri^x 



{1 - fi)qi{x + fi) yU,(x + /i-l) 



3 j2A2{x + /i — 1) 
2 



MbX 



(r2 + r2) 



3/2 



2 (1 - ^^QiV ^2/ 3 ^A2y 

^ y 5 T ^ 

rl To 2 To 



M^y 







(r2 + T^f'^ 
from equations ( [6l [7]) we obtained: 



1/3 



ri = gi 



r2 = 1 + 



'1 - 2rM 1 



3^A2 
2{l-/i) 



2 3 {rl + T2)3/' 

^(1 - 2rjM6 



3 (r2 + T2)3/2 
From above, we obtained: 



g^\2/3 

nV 



1 + 



a; = — /i ± 

;i-2rc)M,(l ^ 



3A2 



3/jA2 



[rl + T2) 



3/2 



-2/3 



X = 1 — /i 



± 



^ ^(1 - 2ro)M, 



-1 -2/3 



-1/ 



1/2 



1/2 



(r2 + r2) 



3/2 
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Figure 2. The figure show the zero velocity curves for /i — 0.025, T = 0.01, in frame one 
qi — 0.00, in second frame show the curves in classical case. 



The triangular equilibrium points are given by putting Q^j. = Qy = 0, y ^ 0, then from 
equations ( [3]) and ( S]) we obtained the triangular equilibrium points as: 



X = —fj, + 



2/3 



'l-2r,)M, 



^ (1-/.) f ^1 ^ 



3 (r2 + T2)=^/2 



2/3 
9i 



4 - + 2 (g^/^ - 2) A, 



4(2r, - l)Mb 



2/3 o 

g/ -3 



2(1^^7^) 



3 (r2 + T2) 



3/2 



1/2 



(12) 



(13) 



All these results are similar with (Szebehely 1967), (Ragos & Zafiropoulos 1995), (Kushvah 
2008 a, Kushvah 20085) and others. 

3. Linearization of the Hamiltonian 



We have to expand the Lagrangian function in power series of x and y, where (x, y) are 
the coordinates of the triangular equilibrium points. We will examine the stability of the 
triangular equilibrium points. For this we will utilize the method of (Whittaker 1965). By 
taking H2-, we will consider linear equations this we have established the relations between 
perturbed basic frequencies. The Lagrangian function of the problem can be written as 



1 n2 



'1 



- /i)gi _^ ii _^ M2 ^ Mb 

n r2 2r| (^^2^2^2)1/2 



(14) 



and the Hamiltonian H = — L + p^x + Pyij, where Px,Py are the momenta coordinates given 
by 

dL dL 

X - ny, py 



Px = jrr 
ox 



= y + nx 

dy 



Let us suppose qi = 1 — e, with |e| << 1 and rewriting the coordinates of triangular 
equilibrium points as: 

2 



x 



2 3 2 3 9 2(1 -/i) 



y 



5e A2 lA^t 4Mb SMfoC /i^zMfoe ^ 
9" ~ 1 9 9 tT ^ 9(1 -/i)/ 

where 7=1 — 2/i. We shift the origin to L4, change x — + x, y 
a = x^ + = y,,,then 



y* 



1 r 2e 2A2e AM^e ^lA^M^ 

— si — — — Ao + r — ~, r 

2\ 3 3 9 (l-/i) 



(15) 

(16) 
y. and 

(17) 



"2~\ ~9"~1 9 9 2^ ^ 9(1 -/i) J 



(18) 



We have to expand L in power series of x and for this we use Taylor's expansion i.e. 

/(x, y) = /(O, 0, ) + [x/,(0, 0) + ?//,(0, 0)] + ^[x2/,,(0, 0) + 2x2//,.,(0, 0) + y V,y(0, 0)] 



+ ^[a:Vxxx(0,0) + 3xVxxy(0,0) + 3x^2/^^^(0,0) +2/3/^j,,(0,0)] + 

Using these values in ( [19]) and with the help of ( [17]) and ( [18]) we get 
L = Lo + Li + L2 + L3 + ■ ■ ■ 
H = Ho + Hi + H2 + H3 + ■ ■ ■ = -L + p^x + pyv 

where Lq, Li, L2, . . . are 

(1 - 
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(20) 
(21) 



Lo= [(a - /i)2 + „2 ^ 



+ 



Va2 + 62 
1 + 



+ 



v/(a-l)2 + 62 ' 2(a-l)2 + 62j ' ((^ _ + 52 + ^2)3/2 

a(l — iJ,)qi 3(a — /i)Mb 



(22) 



Li = X 2n {a — jj) - 

(a — l)/i 
(a - 1)2 + 62)3/2 



^2 + {,2)3/2 _ ^)2 + {,2 + 2-2) 



5/2 



1 + 



3Ao 



26n- 



2 ^(1 - 



2(0- 1)2 + 62 
36Mfe 



a2 + 62)=^/2 _ + ^2 + 2-2) 



5/2 



6^ 



(a -1)2 + 62) 



3/2 



3Ao 



(i;2 + ^2 



2(a- 1)2 + 62 



n 



2 + n{xy — xy) + ^(^^^ + V^) ~ Ex^ — Fy^ — Gxy 



(23) 
(24) 



1728 



2 - 324(2 - 2h^)Mb - 



(2rc - 1) {2 + 15(2 - 7Ai)Mb} Mf, 



(r2 + T^f'^ 
+ 2e {-54 - 270// + 135(10 - 81/x)Mi, 

(2rc - 1) {146 - 33/i + 15(118 - 205//)Mfe} 
(r2 + T2)=^/2 

+ 6A2 {162 - 432/x + 135(2 + 39/i)Mf, 

(2r, - 1) {146 - 240/i - 15(10 - 259/f)Mfe} 



+ 



(r2 + r2) 



3/2 



+ eA^ {144 - 5022// - 270(4 - 395/x)Mb 

(2rc - 1) {458 + 540// + 15(1926 - 7399/t)Mb} 



(25) 



576 



360 + 108(22 + 85//)Mfe + 



(2re - 1) {2 + 5(6 + 35//,)Mb} 



+ 

+ 
+ 
+ 



{rl + T2)3/^ 
2e {54 - 18/( + 45(62 + 309/()M6 
(2rc - 1) {10 + 39/f + 5(202 + lh&h^l)Mb} 

(r2 + T-^fl^ 
6^2 {126 + 45(26 + 99//)Mb 
(2rc - 1) {46 + 5(142 + ^Ql^l)Mb} 

(r2 + 72)3/2 

{576 - 18/t + 90(228 + 1147/t)M6 
(2rc - 1) {522 + 526// + 5(3834 + 27923//)M6} 



+ 



+ 
+ 
+ 
+ 
+ 
+ 
+ 



-1 



(r2 + 72)3/2 



[648 - 1296// + 1620(2 + 3//)Mb 



288^3 

12(2rc - 1) (22 - 44// + 5(34 + 39//)Mb} Mb 

(r2 + T^f^ 
2e {198 - 666/t + 45(94 + 201/t)Mfe 

2(2rc - 1) (190 - 423// + 135(18 + 65//)iV4} iV4 



T2) 



3/2 



6^2 (126 - 468// + 45(26 + 15//)A4 
(2rc - 1) {96 - 260// + 15(66 + 133//)Mfe} 
(r2 + T2)'/' 

6^2 {1368 - 4806// + 90(280 + 429/t)M6 
(2rc - 1) {2106 - 655// + 35(982 + 3025/t)Mb} 



(r2 + T^f^ 



3.1. Perturbed Basic Frequencies 



(26) 



(27) 



Using the method in (Whittaker 1965), to find the canonical transformation from the phase 
space {x,y,Px,Py) into the phase space product of the angle co-ordinates (0i,02) and the 
action momenta co-ordinates Ii, h and of the first order in ll^'^,!^^'^. We consider the 
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following linear equations in the variables x^y: 



-\py 
AX 



dH2 
dy ' 
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The characteristic equation of the Hamiltonian 



Px + p; 



+ n{ypx - xpy) + Ex^ + Fy^ + Gxy 







is given by \A\ = i.e. 

A^ + 2(F + F + n^)\^ + AEF - G'^ + - 2n\E + F) 
Stability is assured only when the discriminant D > 0, where 

D^A{E + F + n^f - 4{4FF - + - 2n^{E + F)} 

from stability condition we obtained: 

/X < /Xco + 0.125885e + {0.571136 + 1.73097e 

0.219964 - 0.398363e + rc(0.202129 + 0.114305e) 

- [0.0627796 - 0.112691e + M^A^ {0.281354 - 1.53665e 
0.654936 - 0.669428e + rc(0.195486 + 0.350878e) 1 

+ irl + r2)^/^ J 

where /Xco = 0.0385209. When D > the roots ±ia;i and ^iu^ {y^x^u^ being the long 
-periodic frequencies) are related to each other as 

ul + [27((1 + n)e - 2) + 9[-18 + 36// + (22 + 69//)] 

54 

12 + 30e + 30// + 95//e) + 135MbA2(18 + 58e + 45// + 188//e) 



+^^-^^^{m, (180 + 2(2re - 1)(44 + 21)//e 
+ 72rc + 4rc^2[78 + 180// + (253 + 873//)e]) } 
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2, ,2 



1^2 



-— [108[e + 3Mb(2 + 7e)] + 18[31e + Mfe(144 + 647e)]A2 
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Mb {4[36 - 47e + Ar^{9 + 28e)] + 3[74 - 297e + 4re(44 + 175e)]A2} 



(r2 + T2) 



3/2 



+ /i 

396(2r, 
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4 8 4 2^M4^ ^ 8 8 



1^ 



4(772rc - 395)e + 12(326rc - 181)^2 + (18452rc - 9667)^26 1 
72 (r2 + T^f^ j 



-27 llle 117 161^26 
+ 



4 8 4 2 
(2r, - 1)[198 - 838e - 1014^2 - 5117^e 



_ ,405 495^2 4185e 25275^26 

Mbi + + + 

^ 4 2 16 16 



36 (r2 + T^f^ 

where Ljj{j = 1,2) satisfy the property (0 < a;2 < < cji < 1). The mass parameter 
is shown in figure < < 2, < ^2 < 1, surface {l)qi = 1 (2) qi = 0.75, 

{3)qi = 0.5,(4)gi = 0.25. 



(34) 



3.2. The Normal Coordinates 

For expressing H2 in a simpler form, we consider the set of hnear equations( [25]) the solution 
of which can be obtained as 

X y Px 

(2nA -G)~ (A2 - n2 + 2E) ~ {nX^ - GX - 2nE + n^) 

^ Py 

(A3 + n^X + 2EX - nG) 

Substituting A = ±iuji and ±za'2, we obtain the solution sets as 

Xj = Kj{2niujj — G), pxj = Kj{—nuj^ — iGuj — 2En + n^) 
Uj = Kj{-ujj -n^ + 2E), pyj = Kj {-iuj^ + iuj{v? + 2E) - Gn] 



(35) 



X 



f+2 = Kj+2{'2niujj - G), Px,j+2 = Kj+2{- 



iGuj - 2En + 



yj+2 = Kj+2[-UJ] 



n^ + 2E), Py, 



J+2 



iLUj{n^ + 2E) - Gn} 



where j = 1,2 and Kj, Kj j^2 are constants of proportionality. Following the method for 
reducing H2 to the normal form, as in (Whittaker 1965), use the transformation 



X = JT where X 
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' Qi ' 
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, T = 


Q2 


Px 




Pi 


. Py . 




P2 



(36) 
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Figure 3. The critical mass ratios in the 3D plot for Mb ~ A^ — [i and q\ = 1, 0.75, 0.5, 0.25, 
rc = 0.9999, T= 0.01 



and J = [Jij]i<i,j<4. i-e. 
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(37) 



Pj = (2/ja;j)^/^ cos 0j, Qi = (^)^''^ sin ^j, (i = 1,2) Under normality conditions: 

XlPx,3 - X3Px,l + yiPyfi - yzPy,l = 1 
X2Px,4 - X4Px,4 + y2Py,A " 2/4Py,2 = 1 

Equivalently, 

- AiuiKiK2{ujl{F -E + n'^) 

+ G^ + + 3^2^ ^ ^2^ _ _ 2n^} = 1 (38) 



- 4iu2K2K4{ujl{F - E + n^) 

+ + 2E'^ + M^E + n^F - 2EF - In"} = 1 (39) 

KjS being arbitrary, we follow the approach of (Breakwell & Pringle 1966) and choose 
Ji,i = Jl,2 



0, implying that 
Ki{2inuJi-G) ■ 



'-{2nuJi — iG), K2{G — 2inu2 



UJ2Ki 



{2nijj2 — iG) 



I.e. 



uji{2nuji — iG) 

K2 



2{2inuji - G) 



hi (say) 
h2 (say) 



U2{2nu)2 — iG) 2(G — 2inu2) 
Using equations ( [3S]),( [Sn]),( SD]) and ( E]) we observe that 

1 



where 



2ujMjMj{M*y 



{u^ -2F + M; = [uj] -2E + r?fl\ 

2{u]-E-2-n^f'\U = ia) 



It is now verified that Ho takes the form: 



H2 
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J=[J, 



ijj4x4 
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(Jl(m|-2n2) ja;2(M|-2n2) 

MiTh M2'M2 
—iUlG iuJiG 



(40) 
(41) 

(42) 



(43) 
(44) 



with I'j 



A/iA/i 
.2 



U12M2M2 
,2 1 U2 



LUlMl 

G 

LuiMiMi 
-nG 

n{uil-Ml) 
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iM2 
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(45) 



4a;| + 9, (j = 1,2) and kl = 2iuf — l.kl = —2ujI + 1. Applying a contact 



transformation from Qi,Q2, Pi, P2 to Q2, P[, P2 defined by (Whittaker 1965) 

a = 1,2) 



dW dW 
Qj 



dP,' 



and 



W 



E 



Q' sin ^ 




I.e. 



cos Pj, 



Pi 



2uj,Q'^smP; (j = l,2) 
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The Hamiltonian H2 is transformed into the form 

H2 = iOiQ[ - UJ2Q2 

Denoting the angular variables P[ and P2 by and 02 and the actions Qi and Q2 by Ji, I2 
respectively, the second order part of the Hamiltonian transformed into the normal form 

H2 = ujih - UJ2I2 (46) 

The general solution of the corresponding equations of motion are 

li = const., (pi = ±uji + const, {i = 1,2) (47) 

If the oscillations about L4 are exactly linear, the equation( |47j) represent the integrals 
of motion and the corresponding orbits are given by 

X = Ji3 A/2a;i/i cos 0i + Jm a/2u^ cos 02 (48) 
y = J2i\ sm 01 + J22 A / sm 02 

\ Ui \ U2 



4. Conclusion 



+ J23 V 2/1 cJi cos 01 + J24 V 2/2^2 sin 02 (49) 



We have seen there are closed curves around the L4(5) so they are stable but the stability 
range reduced due to radiation effect. The effect of oblateness and mass of the belt is 
presented. We have found the normal form of the second order part of the Hamiltonian. For 
this we have solved the aforesaid set of equations. Under the normality conditions, we have 
applied a transformation defined by (Whittaker 1965). We have also utilized the approach 
of (Breakwell & Pringle 1966) for reducing the second order part of the Hamiltonian into the 
normal form. We have found that the second order part H2 of the Hamiltonian is transformed 
into the normal form H2 = Uili — a;2/2- 
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